Magnetic properties of spin-orbital polarons in lightly doped cobaltates 
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We present a numerical treatment of a spin-orbital polaron model for Na^CoC^ at small hole 
concentration (0.7 < x < 1). We demonstrate how the polarons account for the peculiar magnetic 
properties of this layered compound: They explain the large susceptibility; their internal degrees 
of freedom lead both to a negative Curie- Weiss temperature and yet to a ferromagnetic intra-layer 
interaction, thereby resolving a puzzling contradiction between these observations. We make specific 
predictions on the momentum and energy location of excitations resulting from the internal degrees 
of freedom of the polaron, and discuss their impact on spin-wave damping. 
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The complex physics of transition metal oxides of- 
ten leads to intriguing and potentially rewarding prop- 
erties like, e.g., high-Tc superconductivity in the case 
of cuprates or colossal magnetorestistance in the case 
of manganites. The layered cobaltate Na^OoO^ is of 
interest mainly for two reasons: When it is hydrated, 
the resulting NazCoC^yH^O becomes super-conducting 
for doping x ps 0.3, y ps 1.4 1]. Without hydration, it 
shows an exceptionally high thermopower, i.e. the ca- 
pability to transform temperature differences into elec- 
tricity, for 0.5 < x < 0.9 [2], |3J in an unusual combina- 
tion with low resistivity. Magnetic fields strongly affect 
the thermopower in this material Q , which suggests that 
electron-electron correlations and spin degrees of freedom 
are important. These remarkable features of cobaltates 
have triggered large interest due to their potential appli- 
cations, therefore, the origin of the exotic transport and 
magnetic properties of Na^CoO^ is one of the hot topics 
in the field of strongly correlated materials. 

While Na^CoO^ does not show magnetic order either 
for x < 0.7 (except at x = 0.5) || or for x = 1 0, 
it is found in the A-type antiferromagnetic phase below 
T N w 20K for 0.7 < x < 0.9 0, Q. In this phase, the 
carriers move on a hexagonal lattice within the ferromag- 
netic planes, which are antiferromagnetically stacked. 
The lattice spacing within the planes is much smaller 
than between them which typically leads to considerably 
stronger in-plane coupling in such layered systems, e.g. in 
NaNiC>2 with the same crystal structure . Surprisingly, 
this is not the case here: The ferromagnetic in-plane cou- 
pling is only of the order of magnitude of the antiferro- 
magnetic inter-layer coupling [jj or even smaller At 
higher temperatures T ^> Tn, the susceptibility \ shows 
Curie- Weiss behavior, i.e. x — 0), with a negative 

Curie- Weiss temperature 6 < |1 II EE III 03 As has 
been pointed out in Ref. this contradicts the positive 
Curie- Weiss temperature inferred from spin-wave data. 

In this letter, we analyze a spin-orbital polaron model 
proposed in Ref. for NazCoO^ at large x (i.e. small 
hole concentration). By means of unbiased numerical 
techniques, we calculate static and dynamic observables 
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FIG. 1: (Color online) (a) Cubic symmetry for the Co 3+ -ion 
is removed next to a hole (Co 4+ -ion), this reduces the en- 
ergy splitting between t2 g and e 9 levels of Co 3+ ; Hund's rule 
coupling ]h then stabilizes S = 1 instead of S = 0. (b) Vir- 
tual t2 g -e g hopping leading to antiferromagnetic coupling Jij 
between two Co 3+ ions, (c) Spin structure and interactions 
within the polaron: Co 4+ with s — 1/2 in the center, induced 
spin S = 1 on the six adjacent Co 3+ -sites. (Co 3+ -ions further 
away have S = 0.) Parameters: J couples nearest neighbor 
(n.n.) spins with S = 1 on the ring, Jdiag couples them across 
the diagonal, both are antiferromagnetic. J' couples the cen- 
tral s = 1/2 to the ring. 



and thereby provide the first quantitative description 
of the exotic magnetic properties including the negative 
Curie- Weiss temperature 9, the magnitude and temper- 
ature dependence of the susceptibility, as well as a small 
ferromagnetic in-plane coupling mediated by dilute po- 
larons. Further, we discuss the expected impact on mag- 
netic response and predict, apart from the observed spin- 
waves, additional scattering at specific momenta and en- 
ergies, which can be verified experimentally. 

In NaCo02, the Co-d orbitals are split into energeti- 
cally lower t 2g and higher e g states, and the six electrons 
of the Co 3+ -ions occupy the t2 g orbitals, leading to a 
non-magnetic t% g configuration. Consequently, NaCo02 
has rather small susceptibility 0. Removing some of 
the sodium induces holes and the Co 4+ -sites are found 
in a i| g configuration with spin s = 1/2. The few dilute 
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spins embedded in a non-magnetic background, however, 
would give orders of magnitude smaller susceptibility at 
room temperature than observed 0, Hl| . The key to 
the problem is found in the orbital degree of freedom: 
The Co 4+ -holes removes cubic symmetry at the six sur- 
rounding Co 3+ -ions, this splits the £2g-triplet and the e g - 
doublet (on top of the trigonal splitting induced by the 
layered crystal structure) and thus reduces the energy 
gap between the highest t<z g - and the lo wer e g -orbitals, 
see the schematic representation in Fig. 1(a) The gap 



actually becomes smaller than Hund's rule splitting, i.e., 
energy can be gained by transferring one of the electrons 
from the highest i2g-level into the lowest e 9 -orbital and 
forming a triplet with the remaining unpaired electron 

The mechanism of polaron formation by holes induc- 
ing an energy splitting between orbitals is similar to the 
orbital-polaron formation considered for lightly doped 
LaMnC>3 0, however, the 90-degree angle of the Co- 
O-Co bonds causes an important difference: While the 
spin coupling is ferromagnetic in manganites (leading to 
a large total spin of the polaron), the coupling in the 
present case is antaferromagnetic 0|. Because of the 
90° angle, e ff -e g -exchange is suppressed, and the domi- 



nant ti g -e g process, shown in Fig. 1(b) favors antiferro- 
magnetic alignment of the two ions, as well as an analo- 
gous t2 g -t2g hopping. The resulting spin-orbital polaron 
is depicted in Fig. 



1(c) 



A similar antiferromagnetic ex- 
change Jdiag, induced via virtual longer-range hopping 
i', couples the S — 1-sites along the diagonals of the po- 
laron and they also couple to the central s = 1/2. These 
dressed carriers instead of bare holes are in agreement 
with ARPES measurement indicating a strongly renor- 
malized quasi-particle band width |l6l Il7j | . 

The effective spin Hamiltonian resulting from these 
considerations is given by: 



H — J SiSi+l + Jdiag J]] SiSi + 3 + J' SgSi , (1) 



where Si, i = 1, ... ,7, St = Si denote the 5=1 spins 
in the outside ring and sq the s = 1/2 in the center. 
J> ^diagj and J' are the coupling constants: From the 
orbital structure, J « 10- 20meV, Jdi ag < J can be 
inferred 0|. The coupling J' of the central s = 1/2 
to the outside spins is less clear because ferromagnetic 
and antiferromagnetic contributions compete, it can a 
priori be positive or negative. As the Hilbert space of 
this Hamiltonian is only 1458, we easily diagonalize it 
and compute observables at any temperature. 
The susceptibility is particularly interesting: 
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FIG. 2: (Color online) (a) Inverse susce ptibility and (b) sus- 
ceptibility of a polaron, see Fig. |l(c)| for some parameter 
values. For comparison, the susceptibility resulting from a 
single s = 1/2 spin is included in(b) (solid black line without 
symbols) . 



with St ot the z component of the total spin and temper- 
ature T. Ei gives the eigenenergy of state \E{), S? the 
z-component of the spin at site i, fcf, denotes the Boltz- 



mann constant and Z = < 
The inverse susceptibility \~ 



k b T the partition function, 
of a polaron is depicted in 
Fig. 2(a) for several parameter sets. For all cases, the 



slope increases for decreasing T, that is, ((S^ ot ) 2 ) be- 
comes smaller. The reason is the freezing of the S = 1 
system into a singlet state as the antiferromagnetic bonds 
become stronger at T <C J, which leaves only the net 
s = 1/2 spin. If the high-temperature part were contin- 
ued as a straight line, it would cross the T-axis at 8 < 0, 
i.e. give a negative Curie- Weiss temperature. 



i=0 



(2) 



Figure 2(b) shows the susceptibility of a polaron x(T) 
depending on temperature for the same parameter sets. 
We note that comparing the shapes of the curves to ex- 
perimental susceptibility data 0, 0, points towards 
J' < . For comparison, we also included the suscep- 
tibility x = 1/4T of a single s — 1/2, which would be 
obtained, if we had only the s = 1/2 resulting from the 
Co 4+ -site, i.e., if all Co 3+ -ions would remain in the 5 = 
state: For T > J, it is much smaller than the susceptibil- 
ity of the polaron, regardless of parameter values. Exper- 
imentally, the high-temperature values for the suscepti- 
bility Xexp for 0.7 < x < 0.9 range from ~ 10 _4 emu/mol 



FIG. 3: (Color online) Three possible configurations of inter- 
acting polarons. J a t gives the bare ferromagnetic coupli ng be- 
tween s a and sj,, remaining couplings and symbols as Fig.|l(c)] 




FIG. 4: (Color online) Ground-state spin correlations {s a St,} 
between the two Co 4+ sites for the interacting-polaron config- 
urations depending on J', (a), (b) and (c) indicate the cluster 
as depicted in Fig. [3] Jdiag = J, Jab = 0. 



to ~ lCT 3 emu/mol 0, 0, @, EJ El E3 From our high- 
temperature values for the susceptibility of a single po- 
laron Xpoiaron ~ 1/ J", see Fig. |2(b) the susceptibility per 
mole is obtained via 



Xthi 



9^b n a 
J 



(1 X J Xpolaron ? 



(3) 



with Lande g- factor g — 2, Bohr magneton fig, Avo- 
gadro's constant Na and doping x. Taking J = 
20meV, x = 0.82, this corresponds to Xthoory ~ 
1.2 x 10~ 3 , in accordance with the value reported for 
Nao.82Co0 2 Both the shape, see Fig. 2(b) and the 
absolute values of the susceptibility obtained with the 
polaron model therefore agree well with experiment. 

While we have solved the riddle of the susceptibility, we 
have still to show how ferromagnetic interactions within 
the plane can arise and lead to the A-type antiferromag- 
netism observed below 20K To this end, we ana- 

lyze three configurations of interacting polarons depicted 
in Fig. |3 which will occasionally occur as two mobile po- 
larons come near each other. At first, we examine the 



FIG. 5: (Color online) Energies of the lo west states with 
'S'tot = 0,1,2 for the bi-polaron Fig. |3(a)| depending on J' 
for J a t = (main panel) and depending on J a ^ for J' = 
— J(inset). Jdiag = J- 



direct coupling J a b between two Co 4+ -ions situated at 
n.n. sites as in the bi-polaron Fig. 



3(a) 



Like the cou- 
pling between Co 3+ -spins J, it originates from virtual 
hopping processes over a 90-degree Co-O-Co path; it is 
composed of both ferromagnetic and antifcrromagnetic 
terms. With an average hole density n = 1/3 in each t^g 
orbital, we can obtain it from Eqs. (5.5.) to (5.7.) of 
Ref.Hi with A w B « C/2. We arrive at J ab = -8/9A, 
which is ~ — J, because A and J are of the same scale. 
The ferromagnetic net interaction is due to a dominance 
of Hund's-rule coupling over the other terms. 

However, it is not obvious that the more complex po- 
larons do not mediate a different and possibly antifer- 
romagnetic interaction and thus turn the net in-plane 
interaction from ferro- to antiferromagnetic. We there- 
fore solve the complete model given by the clusters in 
Fig. |3 and as the Hilbert space is here larger, we em- 
ploy the Lanczos algorithm. From the low-energy eigen- 
states, we calculate observables like the spin-spin corre- 
lation between the two Co 4+ -sites (s a Sb), which we show 
in Fig.0|in absence of 3 a b- For the bi-polaron Fig. 3(a) 
we see that the ground state favors ferromagnetic align- 
ment, while it is antiferromagnetic for the two clusters 
with larger distance between the Co 4+ -ions. 

In order to make a statement on the size of the result- 
ing effective exchange, we have to consider the involved 
energies. For the bi-polaron Fig. |3(a)[ i.e. for the n.n. 
case, they are shown in Fig. [S] for different values of the 
z-component of the total spin (a good quantum num- 
ber) S* ot — 0,1,2. As expected, the singlet and triplet 
states are degenerate for J' = 0, i.e., when the Co + 
states are not coupled to the Co 3+ system. For \J'\ > 0, 
the ferromagnetic triplet state has lower energy than the 
antifcrromagnetic singlet state, i.e., the effective interac- 
tion via the S = 1-system strengthens the direct ferro- 
magnetic interaction. The effective interaction is then 
combined with the direct coupling J a \, and the inset of 
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FIG. 6: (Color online) Singlet-triplet gap depe nding on J' for 
all three clusters: crosses for the bi-polaron Fig, |3(a)| triangles 
for Fig. |3(b)| and circles for Fig. |3(c)| Jdiag = J 



Fig. shows the energies of the lowest bi-polaron states 
for J' = — J depending on J a t,. We see that the energy 
change is ss Jab(saSb), i-e. the two contributions can be 
simply added J b i po i. « Jbipoi. + Jab- 

Next, we perform a similar analysis on the two other 
configurations of edge- or corner-sharing polarons. As 
can be seen in Fig. 0] we actually find that the ground 
states have antiferromagnetically coupled Co 4+ . The en- 
ergy difference between the lowest triplet state and the 
singlet ground state, which is equivalent to an effective 
antiferromagnetic next nearest neighbor (n.n.n.) cou- 
pling, is depicted in Fig. El While it nearly vanishes 
for J' < 0, it becomes comparable to the ferromagnetic 
n.n. exchange 3 a \, for J' > 0, especially if one consid- 
ers that there are twice as many possible n.n.n. inter- 
actions as n.n. ones (i.e. 12 instead of 6). Such an 
effective antiferromagnetic long-range coupling between 
the polarons should manifest itself as a deviation from 
a n.n. Heisenberg-like dispersion in spin-wave measure- 
ments. While the fact that such a deviation was not 
observed might result from too small momentum vectors 
in experiment, it could also indicate that J' is, in fact, 
negative, which leads to negligible n.n.n. -coupling, see 
Fig. El The spin- wave dispersion thereby provides a way 
to experimentally determine the sign of J'. 

For the n.n. coupling and | J'\ < J, which we believe 
to be a realistic range based on the susceptibility data, 
we read off a ferromagnetic exchange > Jbipoi. ^ —0.5 J 
(see Fig.EJ in addition to the direct J ao ~ — J. However, 
the ferromagnetic coupling is substantially reduced by 
the fact that only some bonds are covered by bi-polarons. 
For hole density S = 1 — x, we have J av ~ (Jbipoi. + 
J a b) ' S. At S Rj 0.2 and J' = — J, this gives an average 
exchange J &v SiSj with J av ~ — 6 meV, which agrees well 
with measured values of ~ — 8 meV |]1LLJ| and ~ — 6 meV 

We have already shown that the internal degrees of 
freedom of the polaron, i.e. the energy-level structure 



of the s = 1/2 coupled to a S — 1 ring with a Hal- 
dane gap, are crucial for the high-temperature behav- 
ior of the static susceptibility x- They likewise influence 
the dynamic susceptibility x"(q,u}), where they lead to 
a quasi-localized mode [l3( in addition to the collective 
and dispersive spin waves. For its experimental obser- 
vation, it is crucial to know energy and momentum of 
this mode. We predict here that its largest weight in 
momentum space is found at q = K = (47r/3, 0,0), a 
little is also obtained for q = M = (0, 27r/v3, 0), but 
the weight is quickly reduced upon approaching the T- 
point (0, 0, 0) (details for the structure factor will be 
given elsewhere) . The energy of the lowest peak is found 
at u> ~ J for Jdiag = J, J' = 0, and is lowered by a 
coupling | J' | > to the central s = 1/2. In spin- wave 
experiments, large damping is expected at these energies 
and momenta because the spin waves strongly scatter 
on the quasi-localized modes. For ferromagnetic J' < 0, 
which has already been noted before to be more plausible 
based on susceptibility and spin-wave data, damping is 
expected to start between uj ~ J/2 and uj ~ J, depending 
on J', as actually observed in neutron-scattering experi- 
ments []j,|8j]. For an antiferromagnetic J' > 0, the energy 
of the localized mode is reduced much faster and reaches 
uj = for J' = J. This seems to rule out large anti- 
ferromagnetic J', because spin- waves are actually clearly 
measurable at not too large energies 0, • 

To summarize, we present a coherent and quantita- 
tive description of the magnetic properties of Na^OoC^ 
at large x. The concept of spin-orbital polarons could 
likewise be a key in understanding the unconventional 
transport properties of cobalt-based materials. 

We would like to thank B. Keimer, C. Bernhard, S. 
Bayrakci and A. T. Boothroyd for useful communica- 
tions. 
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